A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

12. (08.14.1; 08.23.1; 02.04.1; 02.03.2) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Transport properties of degenerate electrons in neutron star 
envelopes and white dwarf cores 



A.Y. Potekhin 1 *, D.A. Baiko 1 , P. Haensel 2 , D.G. Yakovlev 1 

1 Ioffe Physical- Technical Institute, Politekhnicheskaya 26, 194021 St .-Petersburg, Russia 

2 N. Copernicus Astronomical Center, Bartycka 18, PL 00-716 Warsaw, Poland 

Received 3 December 1998 / Accepted 2 March 1999 



in 

m 
> 

(N 

m 
O 

On 

Or 

6 
to 

C3 



Abstract. New calculations of the thermal and electrical 
electron conductivities are performed for a broad range 
of physical parameters typical for envelopes of neutron 
stars and cores of white dwarfs. We consider stellar mat- 
ter composed of astrophysically important chemical el- 
ements from H to Fe in the density range from 10 2 - 
10 4 gem -3 up to 10 7 -10 10 gem -3 , where atoms are fully 
ionized and electrons are strongly degenerate. We have 
used modified ion structure factors suggested recently by 
Baiko et al. (1998). In the ion liquid, these modifications 
take into account, in an approximate way, instantaneous 
electron-band structures that reduce the electron-ion scat- 
tering rate. In crystallized matter, the new structure fac- 
tors include multi-phonon processes important at temper- 
atures not very much lower than the melting temperature 
T m . The transport coefficients obtained differ significantly 
from those derived earlier in the important temperature 
range T m /5 < 5T m . The results of our numerical cal- 
culations are fitted by analytical expressions convenient 
for astrophysical applications. 
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1. Introduction 

Thermal and electrical conduction in the envelopes of neu- 
tron stars and the cores of white dwarfs plays crucial role 
in many aspects of evolution of these stars. Thermal con- 
ductivity is the basic quantity needed for calculating the 
relationship between the internal temperature of a neu- 
tron star and its effective surface temperature; this rela- 
tionship affects thermal evolution of the neutron star and 
its radiation spectra (e.g., Gudmundsson et al. 1985 ; Page 



1997; Potekhin et al. 1997, hereafter Paper I). Electrical 



conductivity is the basic quantity used in calculations of 
magnetic-field evolution in neutron star crusts (e.g., Mus- 
limov & Page 1996; Urpin & Konenkov 1997; Konar & 



Bhattacharya 1997). Thermal conductivity of degenerate 
matter is also an essential ingredient for white-dwarf pul- 
sation modelling (e.g., Fontaine & Brassard 1994). 

For applications one should know the transport prop- 
erties of dense stellar matter where electrons are strongly 
degenerate and form nearly ideal Fermi-gas, whereas ions 
are partially or fully ionized and form either a strongly 
coupled Coulomb liquid or a Coulomb crystal. Under such 
conditions, electrons are usually most important charge 
and heat carriers, and the electrical and thermal conduc- 
tivities are mainly determined by electron scattering off 
ions (hereafter, ei scattering). 

The conductivities of degenerate electrons due to ei 
scattering were studied in a number of papers. The general 
formalism, based on a variational metho d (Zi man 1960), 
has been developed by Flowers & Itoh ( 1976|) (see refer- 
ences to earlier results therein). Their numerical results, 
however, have been critically revised by Yakovlev & Urpin 
( 1980| ), who developed a simple analytical description of 
the conduction due to ei scattering in dense Coulomb plas- 
mas. The results by Yakovlev & Urpin ( |1980 ) were con- 
firmed in detailed calculations for both solid (Raikh & 
Yakovlev 1982) and liquid (Itoh et al. 1983| ; Nandkumar 
& Pethick p"984| ) regimes. Later Yakovlev ( |1987| ) calcu- 
lated the electron transport coefficients in the liquid phase 
taking into account non-Born corrections. In Paper I, the 
authors performed extensive calculations of the same co- 
efficients including additionally the effects of responsive 
electron background on the ion structure factors. 

Itoh et al. (1984, 1993) improved the results by 
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Yakovlev & Urpin fll980| ) and Raikh & Yakovlev ( |1982[ ) in 
the solid regime by including the nuclear form factor F(q) 
(where hq is a momentum transfer), in order to take into 
account finite sizes of atomic nuclei, and studied the role of 
the Debye-Waller factor e~ 2W ^ (e.g., Kittel |l963| ) , which 
describes reduction of ei scattering rate in a crystal due to 
lattice vibrations. The Debye-Waller factor proved to be 
important at temperatures close to the melting temper- 
ature of a Coulomb crystal, and at sufficiently high den- 
sities, where zero-point vibrations are large. Later Baiko 
& Yakovlev ( |l995[ |l996[ ) made detailed Monte Carlo and 
analytical calculations of the electron transport in crys- 
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where n e is the electron number density, m c is the electron 
mass, x r w 1.009 (peZ/A) 1 / 3 is the relativistic parameter, 
Z is the ion charge number, A is the atomic weight, and 
Pq is mass density p in units of 10 6 gcm~ 3 . The electron 
degeneracy temperature is 



(e F - m c c 2 )/fc B w 5.93 x 10 9 (V 1 + ^ 2 - l) K, (2) 



where /c 



B 



IS 



the Boltzmann constant and ep = 
x% = m*c 2 is the electron Fermi energy. Our 
analysis is thus limited by the condition T <C Tp. It is also 
restricted to temperatures T^5x 10 9 K at which atomic 
nuclei do not dissociate. Electrostatic screening properties 
of the electron gas are characterized by the Thomas-Fermi 
wave number fcTF: 



fc TF = 47re 2 



dn c 
dp 



Off y/1 + X\ 



(2k F ) 2 



where p ~ ep is the electron chemical potential and ctf = 
e 2 /he = 1/137.036 is the fine-structure constant. 

The state of ions (atomic nuclei) can be conveniently 
specified by the Coulomb plasma parameter 



r = 



1/3 



talline dense matter including the nuclear form factor and 
the Debye- Waller factor; they fitted the results by simple 
analytical expressions. These results were in reasonable 
agreement with those obtained by Itoh et al. (1984, 1993[ ). 

Nevertheless the transport theory developed in the 
cited articles possessed one important drawback: it pre- 
dicted unrealistically large jumps (by a factor of 2-4) of 
the transport coefficients at the melting point in spite of 
the fact that many other physical properties of ion liquid 
and solid were very similar. This indicated that the theory 
was incomplete and had to be improved. 

The improvements have been suggested recently by 
Baiko et al. ( 1998j ) (hereafter Paper II). In the solid 
regime, multi-phonon processes have been included in the 
electron-phonon scattering, whereas all previous calcula- 
tions used the onc-phonon approximation. Furthermore, it 
has been noted in Paper II, that the static structure fac- 
tor of ions conventionally used in the liquid regime may 
require modification when applied to calculation of the 
electron scattering processes because of appearance of in- 
cipient ordered structures. The authors suggested an ap- 
proximate treatment of this effect by subtracting a certain 
part from the static structure factor. Both modifications 
affect significantly the electron transport properties near 
the melting point and reduce the jumps of the transport 
coefficients. 

In this paper, we apply the results by Baiko et al. 
( |1998| ) to calculation of the electron electrical and ther- 
mal conductivities in a wide range of physical parameters 
typical for the envelope of a neutron star or the core of a 
white dwarf. We also derive an effective ei scattering po- 
tential that yields the conductivities in a simple analytical 
form, reproducing the numerical results within, at most, 
a few tens of percent in the whole of the physically mean- 
ingful range of the plasma parameters. Finally, we discuss 
the main features of the electron transport properties, and 
the role of various electron scattering mechanisms. 

The paper is organized as follows. In Sect. || we in- 
troduce basic definitions and give a brief overview of the 
main features of electron conduction in different regimes. 
In Sect. H we describe our calculation of the electron trans- 
port coefficients and propose a fit to these coefficients. 
Numerical results are discussed in Sect. |[ 

2. Dense degenerate matter 

2.1. Equilibrium properties 

Consider fully ionized degenerate stellar matter in the den- 2 - 2 - Transport coefficients and structure factors 
sity range from about 10 2 -10 4 gem -3 to 



(3) 



(Ze) 2 _ 22.75 Z 2 
fc B Tai T 6 



(4) 



where e is the elementary charge, ai = [3/ (AnnC)] 1 / 3 is the 
ion-sphere radius, — n c /Z is the number density of ions, 
and Tq is temperature in units of 10 6 K. If T <C 1, ions are 
weakly coupled and form the Boltzmann gas. For T > 1, 
they constitute a strongly coupled liquid. Freezing occurs 
at a temperature T m which corresponds to T = T m . For 
classical ions, F m = 172, whereas quantum effects (zero- 
point ion vibrations) suppress the freezing and increase 
r m (Nagara et al. 1987). The quantum effects become im- 
portant at T <C T p , where 



T p = hu p /k B w 7.832 x 10 6 (Z/A)pl /2 K 



(5) 



is the ion plasma temperature, u p = (AnZ e ni/mi) ' is 
the ion plasma frequency, m\ — Am u is the ion mass, and 
m u = 1.66 05 x 10 -24 g is the atomic mass unit. Under 
realistic conditions, the quantum effects strongly suppress 
crystallization of hydrogen and helium plasmas, but do 
not affect significantly the melting of carbon and heavier 
elements (e.g., Chabrier 1993 ). However, they affect the 
properties of matter of any composition at T <C T p . 



10 10 gcm~ 3 . 

For simplicity, we assume that there is one ion species at 
any given p and T. 

The state of degenerate electrons can be described by 
the Fermi momentum p-p or wave number kp: 



Electrical (a) and thermal (k) conductivities of degenerate 
electrons can be conveniently expressed through effective 
electron collision frequencies, v a and i> K , as (e.g., Ziman 
1960|; Yakovlev & Urpin |1980D 



PF 



Hkv, 



h (37r 2 n e ) 1 / 3 = m e cx T , 



(1) 



n c e" 
m*v„ 



7T 2 fc|Tn c 
3to*i/ k 



(6) 
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The collision frequencies are reduced to sums of partial 
collision frequencies associated with relevant electron scat- 
tering mechanisms which can be studied separately. This 
approximation is accurate to ~ 1% in the case of strongly 



degenerate electrons (e.g., Ziman 1960; Lampe 1968). 
In the solid phase at T < T p , where the frequency of 



collisions 



is strongly reduced, the electron trans- 



port is limited by scattering off various irregularities of the 
crystalline structure. The particular case of ion impurities, 
occasionally embedd ed in t he lattice, was studied in detail 
by Itoh & Kohyama (|l993|) and will be taken into account 
in this article; in this case charge and heat transport are 
determined by a single collisional frequency fi mp . On the 
other hand, in the liquid phase, at not very strong elec- 
tron degeneracy for Z ^ 6 the thermal (but not the electri- 
cal) conductivity may be affected by electron-electron (ee) 
collisions. The ee collision frequency v ee was evaluated, 
e.g., by Urpin & Yakovlev ( |l980[ ), Timmes (|l992| ), and 
in Paper I [see Eq. ([52]) below] . These results for electron- 
electron and electron-impurity scattering are not modified 
by the present consideration. The total effective collision 
frequencies are v a = v e J , i/ K = j/f + vff in the liquid and 
v a , K = v™ K + u imp in the solid. 

We will focus on ei scattering. In a weakly coupled 
ion gas, T <C 1, collective effects lead to a dynamical ion 
screening of ei interaction which can be described by the 
dynamical dielectric function formalism (e.g., Williams & 
DcWitt 1969). In a strongly coupled ion liquid, it is cus- 



tomary to use the static structure factor of ions in order 
to describe the correlation effects (Hubbard 1966] ). This 
description, however, does not apply to quantum liquids, 
such as H or He at high densities (e.g., Chabrier 1993| ). In 
crystalline matter, ei interaction is adequately described 
in ter ms of absorption and emission of phonons (Abrikosov 



1961 ). The description can be realized using a dynamical 



structure factor of ions (Flowers & Itoh 1976 ). 

The ei collision frequencies can be expressed through 
dimcn sionlc ss Coulomb logarithms A ffiK (cf. Yakovlev & 
Urpin [19801) : 



plv F 



•A, 



4Ze F 
3irh 



af A CTiK , 



(7) 



where v F = Pf /w* is the electron Fermi velocity. 

For a strongly coupled plasma of ions (r <; 1), the 
Coulomb logarithms calculated in the variational ap- 
proach (with the simplest trial functions, Ziman 1960| ) in 
the Born approximation read 



A„ 



2I, V 



dqq u (q) S a . K (q) 



2k F J 



(8) 



where qo is the cutoff parameter, equal to zero for the 
liquid phase and to the equivalent radius of the Brillouin 
zone <7b = (6-7r 2 n i ) 1 / 3 in the solid phase. The latter cutoff 
filters out umklapp electron-phonon processes, which op- 
erate at q > qB and give the main contribution to A CT K , 



from normal processes, that take place at q < q^ and 
are negligible under the conditions of study (e.g., Baiko 
& Yakovlev p995| ). Furthermore, u(q) = \U(q)\/{4irZe 2 ), 
U{q) is the Fourier transform of the elementary ei scat- 
tering potential, the factor in square brackets describes 
kinematic suppression of backward scattering of relativis- 
tic electrons (e.g., Berestetskii et al. 1982), and S G ^{q) 



are the effective static structure factors which take into 
account ion correlations, as discussed below. 

In order to take into account corrections to the Born 
approximation, we multiply additionally the integrand in 
Eq. (||) by the ratio of the exact cross section of Coulomb 
electron scattering to the Born cross section. This approx- 
imate treatment of non-Born corrections was proposed by 
Yakovlev ( 1987 ) and used in Paper I. The corrections are 
significant for Z J> 20 and p i> 10 6 gcm~ 3 . 

In the case of Coulomb scattering, one has 



u(q) = 



F(q) 
f\e(q)\' 



(9) 



where e(q) is a static longitudinal dielectric function of the 
electron gas (Jancovici 1962| ), describing electron screen- 
ing of the scattering potential. The static (ui — > 0) approx- 
imation of electron screening is adequate as long as typi- 
cal energies transferred hu are small compared with hqvp , 
which is the case, since the momentum transfer fiq ~ p-p. 
At densities p < 10 10 g cm -3 and under the condition of 
full ionization, one can safely set the ion form factor F(q) 
equal to 1, which corresponds to point-like nuclei. 

The structure factors in Eq. (M) are given by (Paper II) 



Sa{q) 

S K (q) 
SS K (q) 



dcuS"(q,Lu) 



S<j(q) 



3fcp 



8S K (q), 



dcuS"(q,uj) 



1 



(10) 

(11) 
(12) 



In this case z — TiLo/lk-sT) and S"(q,ui) is the inelas- 
tic part of the total dynamical structure factor S(q,Lu) = 
S"(q,iv) + S'(q)5(u), whereas the elastic part, S'(q)S(uj), 
describes Bragg diffraction. When interaction of electrons 
with a crystalline lattice is considered, the Bragg diffrac- 
tion leads to appearance of electron band structure (Bloch 
states) but does not c ontribute to electron transport (e.g., 
Flowers & Ito h h976|) . If z values "allowed" by S"(q,oj) 
in Eqs. (|Io|)-(|l2) are small, as it happens for scattering 
in a classical Coulomb system [i.e., at T <; T p ; cf. Eq. 
(|l8|)1, we can pull the other functions containing z out 
of the integral setting z — 0. Then SS K (q) vanishes, and 
S a (q) = S K (q) = S"(q), where S"(q) = S"(q, u) 
is the inelastic part of the static structure factor. In this 
case the variational solution that is used in Eq. (^|) be- 
comes exact. 

In the liquid phase, T < r m , it is only the static struc- 
ture factor in the classical regime S(q) = f^°° duiS(q, lo) 
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(e.g., Young et al. 1991 and references therein) that has 
been determined quite accurately. Thus we will consider 
only classical liquids, T <; T p . In the solid phase S(q,u>) 
was calculated with reasonable accuracy in Paper II, which 
enables us to study the transport properties of quantum 
and classical solids. 

3. Calculation of electron conductivities 

3.1. Solid phase 

The dynamical structure factor of a Coulomb crystal has 
been calculated in Paper II in the harmonic-lattice ap- 
proximation, taking into account explicitly multi-phonon 
processes. For q > the inelastic part may be written as 



S"(q,u>) 
K(q,T,t) 



-2W{q)-hul/(2k B T) 



dte~ iuJt 



exp 



2tt 
kg 2 

2m\ \ io v sinh (z v /2) 



COS UJyt 



K(q,T,t),(13) 
1, (14) 



where (. . .) p h denotes averaging over the phonon spectrum 
in the first Brillouin zone, and 



W(q) = 

2771; 



(15) 



P h 



In this case, v = (Q,s), s = 1,2,3 enumerates phonon 
polarizations, Q is a phonon wave vector, e„ the polar- 
ization vector, u v the frequency, and n v = (e z " — 1) 1 is 
the mean number of phonons, z v = Tiu) v / '(ksT). For the 
lattice types of interest [e.g., body centred cubic (bcc) or 
face centred cubic (fee) ones], W(q) = r^q 2 /6, where r\ is 
the mean-squared ion displacement. Thus W(q) does not 
depend on the orientation of q. An analytical fit to W(q) 
was proposed by Baiko & Yakovlev ( 1995| ): 



W(q) = -a, 



U-\ e 



-9.177 



2rj 



(16) 



where 77 = T/T p , u n = ((uv/w p ) n ) p h is a frequency mo- 
ment of Coulomb lattice (tt-2 = 13.0 and u_i = 2.8 for 



bcc lattice, cf. Pollock & Hansen 1973), and 



OL\ = Q() 



a 



3Tr] 



1.683,/-^. (17) 



It is possible to derive an asymptote of S"(q,u>), 
Eq. @, for T > T p (classical solid) and for \tu\ > u p . 
Expanding cosuj^t «1- (uj l/ t) 2 /2 and sinh {z v /2) « z v j2 
and noting that the second term in Eq. (|lj) does not con- 
tribute into S"(q,w) at these frequencies, we obtain 



S"(q,u) 



l"K qvT 



■ exp 



q 2 v\ 



Tito 
2fc B T 



(18) 



where vt — \/2k^,T ~Jm\ is the thermal ion velocity. Thus 
in a classical solid the collisional energy transfer is limited 
either by TiqvT or by fkj p , in both cases being smaller than 



k^T. Therefore only the values z 1 contribute to the 
integrals (|i"o|), ( |l2| ) in this case. 

In the general case of classical or quantum Coulomb 
crystals, the effective structure factors (|l0|) and (|l2|) can 
be written as (Paper II) 



n t 

cosn x 



(19) 



5S K (q) 



-2W(q) 



Ax 



1 — 2 sinh x 
cosh 4 x 



K(q,T,t), (20) 



where x — irtT/Ti. We have calculated S a (q) and SS K (q) 
for the bcc lattice and fitted them by the expressions 



SS K (q) 
Wi(q) 



e -2W(q) ^ e 2Wi( 9 ) _ ^ 



Oil 



91?, 2 e 



-2W(q) 



(1 + lllAr) 2 ) 2 

0.101 T] 4 



(0.06408 + 7y 2 )(0.001377 + iff/ 2 

o2 



(21) 



(22) 



2y/ V 2 + (U- 2 /117) 2 



These fits cover a wide range of the parameters, 0.001 < 
77 < 10 and < «i < 0.3, sufficient for calculation of 
transport coefficients. Equations (21) and ( |22| ) reproduce 
also the asymptotes of the effective structure factors at 
low and high rj which can be obtained from Eqs. (|l^) and 
(^o|). The maximum fit error, equal to 4%, occurs for ct\ — 
0.001, and r] = 0.04. 

We have also calculated the effective structure factors 
for fee Coulomb lattice and the results appear to be al- 
most indistinguishable from those obtained for bcc lattice. 
Therefore, the electron transport coefficients are insensi- 
tive to the lattice typeQ and we will calculate them for 
bcc lattice. 

3.2. Liquid phase 

In the classical liquid, we employ the static structure fac- 
tor, S(q), obtained by Rogers and DeWitt (unpublished) 
for the one-component classical plasma of ions in a rigid 
electron background by solving the modified hypernetted- 
chain (MHNC) equations (Ros enfeld & Ashcroft |1979|) , 
and fitted by Young et al. Jl99l[ ) in the range 1 < T < 225. 
For calculating ei scattering rates, we have modified these 
structure factors by subtracting the contribution of elastic 
scattering as prescribed in Paper II: 



S<r, K (q) = S{q) 



G#0 



(23) 



1 The results by Baiko & Yakovlev (1995, 1996) for the fee 
lattice, which led to a different conclusion, were inaccurate 
due to an error in a Brillouin zone integration scheme for this 
lattice. 
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where the sum is over all non-zero reciprocal lattice vec- 
tors G, and the upper bar means averaging over orien- 
tations of q. This modification is meant to account for 
instantaneous electron band-structures which emerge in a 
strongly coupled Coulomb liquid because of local tempo- 
rary crystal-like ordering. In this context, the choice of 
the lattice type and corresponding vectors G is ambigu- 
ous, but the Coulomb logarithms are insensitive to it. 

Since our consideration is based on the classical static 
structure factor S(q), we obtain A CT = A K and v e a l = 
in the liquid regime. This is justified because we typically 
have T ^ T p , for the cases under study (see below). Nev- 
ertheless we could easily incorporate the quantum effects 
into the calculation, were the quantum dynamical struc- 
ture factors S(q, u) available for ion fluid. 

Note that in Paper I the Coulomb logarithms in the ion 
liquid were calculated with S(q) obtained for a polarizable 
electron background including the local field corrections. 
The effect of the response of the background appeared to 
be noticeable for H and He plasmas only. We neglect this 
effect in the present calculations because our consideration 
of ion liquid for light elements is approximate anyway due 
to the neglect of quantum corrections to S(q). 

For light elements, the highest temperatures corre- 
sponding to the liquid regime, T > 1, are much below 
T-p. Accordingly, there exists a temperature range where 
T < Tp and T < 1 . In this interval the formalism of the 
effective structure factors does not provide an accurate 
treatment of ion screening. For T < 0.25, the Coulomb log- 
arithms were calculated in Paper I taking into account dy- 
namical character of ion sc reeni ng at T <C 1 and fcxF <C kp 
(e.g., Williams & DeWitt |l969| ): 



A„ 



?D / 2 
^|/(2c 2 )+/3/2, 



1 + C 



In 



1 



i + C 



(24) 



qr> = vST/di is the inverse Debye screening length of ion 
gas,0 C = (^tf/<Zd) 2 , and (3 = naiZvp/c; (3/2 is a lowest- 
order non-Born correction in the weak electron-screening 
approximation (Yakovlev 1987 ). 

The Coulomb logarithms in the transition domain from 
weak (r <C 1) to strong (r ^ 1) ion coupling can be cal- 
culated using the formalism by Boerker et al. (1982). We 
do not apply this formalism, but the Coulomb logarithms 
calculated at T < 0.25 and at T > 1 converge nicely and 
can be fitted in a unified manner. This convergence de- 
teriorates at lower p ~ 10-100 gem -3 , because electron 
screening ceases to be weak (/ctf ~ ^f), and Eq. ( p4|) 
becomes inaccurate. 



2 Definitions of the screening parameters below Eq. (15) of 
Paper I, corresponding to the present Eq. (p4|), contained sev- 
eral misprints. Correct definitions are reproduced here. 



3.3. Numerical results and fitting formulae 

Using the effective structure factors S a ^(q) described 
above, we have calculated the Coulomb logarithms A CTiK 
for Z from 1 to 26 and for mass numbers A corresponding 
to the most abundant isotopes. The mass density p var- 
ied from 10 gem -3 for Z = 1,2 and from 100 gem -3 for 
Z > 3 to 10 10 gem -3 ; the coupling parameter T varied 
from 1 to 10 4 for Z < 20 and to 10 5 for 20 < Z < 26. The 
physically meaningful domain of the parameters is con- 
strained by several conditions. First, it is assumed that 
the atoms are fully ionized (for a treatment of the case of 
partial ionization, sec below). Secondly, light elements are 
not present at very high densities and temperatures since 
they burn into heavier ones. Thirdly, our calculation in the 
liquid phase is confined to the classical regime (T > T p ), 
where the quantum corrections to the ion structure factors 
are neglected. Fourthly, electrons are assumed to be degen- 
erate (T <§; Tp). Finally, the present formalism appears to 
be invalid at very low temperatures, T <C T p Z 1 / 3 e 2 / (hvp), 
where the electron band-structure effects strong ly reduce 
ei scattering rate (e.g., Raikh & Yakovlev 1982), which is 
not taken into account in Eq. (||). 

For practical applications, especially for modelling 
thermal and magnetic evolution of the neutron stars, as 
well as pulsations of the white dwarfs, it is desirable to 
have analytical formulae for the transport coefficients, 
rather than tables, graphs or cumbersome theoretical ex- 
pressions. Analytical formulae for a and n were presented 
in several papers based on the earlier theoretical results 
described in Sect. (Flowers & Itoh tL98lt Yakovlev & 



Urpin |1980[ Itoh et al. |1993| ; Baiko & Yakovlev [1995| Pa 



per I). As our present results are significantly different, we 
propose new analytical expressions, which combine rea- 
sonable accuracy with simplicity. 

Instead of constructing ad hoc fits to the numerical 
values of a and k, we have chosen to devise an effective ei- 
scattering potential that would allow us to perform explicit 
analytical integration in Eq. (^J) and that would reproduce 
correctly the familiar limiting cases: the case of Debye 
ion screening in a weakly-coupled plasma and the case of 
scattering by high-temperature phonons. In the first case, 
u 2 (q)S(q) in Eq. (||) should be replaced by (q 2 + <? 2 ) -2 , 
where q s is the inverse screening length. In the second 
case (at T p ^ T < T m ), ei scattering rate is determined 
by u 2 {q)S"{q), where u 2 {q) ra (q 2 + fc 2 F ) -2 , and S"(q) is 
the approximate effective static structure factor (Paper II) 
which can be written as S"(q) « 1 — exp[— u-ia 2 q 2 j (3r)]. 
This approach ensures that the analytical limits men- 
tioned above are reproduced automatically by the fit ex- 
pression. 

We propose the following form of the effective potential 
in Eq. (§: 



[u 2 (q)S a , K {q)] eS 



1 



(q 2 + q 2 ) 2 



(25) 
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where e~ w ^ plays role of an effective Debye-Waller fac- 
tor at large T and is negligible at T < 1, q s is an effective 
screening wave number, D{rf) is associated with the quan- 
tum correction to the Debye-Waller factor, and G a>K {r), 0) 
is a phenomenological factor that describes reduction of 
thermal ion displacements in quantum solid at T ^ T p 
and contains non-Born corrections expressed through the 
argument [see Eq. (p4|)] . Our numerical results are repro- 
duced with the following choice of q s , w(q) and G atK (r], 0): 



r 



<z. 2 (r) 

w(q) 



= ^(l + 0.06r)e- vT , 
= U - 2 (g/<ZD) 2 (l + /3/3), 



(1 + 0.122/3 2 ), m 



0.19 



^fW Q 1 - ZV6' 

G a (r},0) + 0.0105 (l-Z- 1 ) 



V 



D{rj) = cxp(- 



otQU-i e 



(rf + 0.0081) 3 / 2 : 
9.m /4 ) ; 



where ao is given by Eq. (|17|). Inserting Eq. ( |25| ) into 
Eq. (||), we arrive at the analytical expressions for the 
Coulomb logarithms presented in the Appendix. 

The typical error of our fits in the physically reason- 
able range of parameters is 3% (maximum 6%) for Z > 20 
and gradually increases with decreasing Z . The maximum 
error occurs for low Z at the melting point in the high- 
density region, where our formulae interpolate across the 
conductivity jump discussed in the next section. For in- 
stance, a typical error for carbon plasma is 8% and the 
maximum is 22% at the highest density (p = 10 9 gem -3 ) 

and r = r m . 

Finally, let us outline the case of multicomponcnt ion 
mixtures. Actually the case deserves a separate study and 
we treat it approximately here. At least for T <; T p , it 
would be a good approximation to replace Z 2 riiA e i — > 
Y]j Z 2 njA e j in Eq. (j?j), where summation is over all ion 
species j, and the Coulomb logarithm A e j depends gen- 
erally on j. In a strongly coupled ion system we recom 
mend to calculate A e j 



z 



Zj and r,- 



from Eqs. flAT] ), @-@ with 
(4 7 rn /3) 1 / 3 (fc B T)- 1 (the ion- 
coupling parameter for ion species j). The latter expres- 
sion is prompted by the "additivity rule" that is highly 
accurate for thermodynamic functions of multicomponcnt 



ion mixtures (Hansen et al. 1977) 



Another option is to adopt the widely used mean-ion 
approximation. In the latter approximation, the plasma 
is treated as a mixture of electrons and one ionic species, 
with an effective charge eZ e g equal to an average charge 
of all ions at different ionization stages. The mean-ion ap- 
proximation can be used also in the regime of partial ion- 
ization (cf. Paper I). 
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Fig. 1. Dependence of the electrical conductivity pro- 
duced by ei scattering on temperature (lower horizontal 
scale) or Coulomb coupling parameter (upper horizontal 
scale) in carbon plasma at p = 10 4 g cm~ 3 . Filled circles 
show the present numerical results, and the solid line is 
given by our fitting formula. The dashed line is obtained 
under traditional assumptions (see the text) and exhibits 
a jump at the melting point. Dotted lines are obtained 
including electron scattering by 16 O impurities with con- 
centrations xi mp = 0.02 and 0.04. 



4. Discussion of the results 

Fig. |l| shows the temperature dependence of the electrical 
conductivity of degenerate electrons in a carbon plasma 
at p — 10 4 g cm -3 . Fig. || shows the same dependence 
of the thermal conductivity. Upper horizontal scale indi- 
cates corresponding values of the ion coupling parameter 
F. Since the ion charge number is rather low, Z = 6, non- 
Born corrections are insignificant. All the data presented 
in the figures, except dotted lines, show the conductivities 
produced solely by ei scattering. 

Filled circles display our present numerical values of 
the conductivities, while solid lines are the analytical fits. 
Dashed lines show the 'old' conductivities calculated us- 
ing the approximations which have been widely employed 
in the previous works (Sect. |l|): the one-phonon approx- 
imation in the solid phase and use of the total (inelastic 
+ elastic) ion structure factor in the liquid phase. One 
can see large jumps of the 'old' curves at the melting 
point. These jumps were present for all elements and for 
all plasma parameters, and they were typically a factor of 
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Fig. 2. Same as in Fig. [I] but for the thermal conductiv- 
ity. 



2-4 in magnitude (Itoh et al. 1993). The modification of 
the structure factors improves the treatment of the con- 
ductivities both in the liquid an d soli d ph ases of strongly- 
coupled ion system (Sects. 3T and |3.2| ) and makes the 
jumps almost invisible for different chemical elements in a 
broad range of densities. This has allowed us to produce 
the unified fits which are equally valid in solid and liquid 
matter. 

Nevertheless, our calculations do show large jumps of 
the conductivities at the melting point at high densities, 
where zero point vibrations become important. We sug- 
gest that these jumps are artificial and come from using 



classical ion structure factors in ion liquid (Sect. 3.2) un- 



der the conditions in which quantum effects in liquid are 
really important. On the other hand, the quantum effects 
are properly included in our calculations for crystalline 
matter. Since the numerical data used for constructing 
the fitting formulae include both phases, liquid and solid, 
the fitting of these data by the unified analytical expres- 
sions shifts the conductivities in the liquid phase closer to 
those in the solid phase. Therefore we expect that the fits 
in the high-density ion liquid give more reliable electron 
conductivities than our original numerical data. This as- 
sumption will be checked in the future when the structure 
factors in ion liquid are calculated taking into account the 
quantum effects. 

Dotted lines in Figs. |l| and ^Jshow the effect of another 
electron scatter ing mechanism - scattering by charged im- 
purities (Sect. 1^2]). We have assumed an admixture of 



Fig. 3. Electrical conductivity produced by ei scatter- 
ing in Fe matter at p = fO 8 g cm~ 3 vs temperature 
or Coulomb coupling parameter. Filled circles show the 
present numerical results, and the solid line is the fit (non- 
Born corrections included). Open circles are the present 
results but in the Born approximation. The dashed line is 
the old result in the Born approximation. 



oxygen nuclei with concentrations 2% and 4%. Electron- 
impurity scattering is seen to have little effect on the con- 
ductivities at high temperatures, but it becomes dominant 
scattering mechanism at T <^ T p . 

Figs. H and || are analogous to Figs. |f| and ^ and 
show the temperature dependence of the conductivities 
produced by ei scattering in iron plasma at p = 10 8 g 
cm~ 3 . Filled dots are our numerical results and solid lines 
are the fits. One can again see large jumps of the tradi- 
tional conductivities at the melting point, and the smooth 
character of the improved curves. For elements with high 
Z (like Fe), the non-Born corrections become important 
in a dense plasma (Yakovlev 1987). To illustrate this ef- 
fect, open circles display results of our calculations in the 
Born approximation. The non-Born corrections are seen to 
increase the ei collision frequency (decrease the conduc- 
tivities) by about 20-30%. The corrections become lower 



when density decreases below 10 g cm . Dashed lines 
show the 'old' conductivities calculated in the Born ap- 
proximation. The divergence between the 'old' and new 
results in the Born approximation, seen in Fig. ^ at rela- 
tively low temperatures, is caused by a not very accurate 
determination of the low-temperature (T <C T p ) asymp- 
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Fig. 4. Same as in Fig. 
ity. 



but for the thermal conductiv- 



tote of an effective electron-p honon potential in the 'old' 
results by Baiko & Yakovlev ( 1995 ). 

The density dependence of the thermal conductivity 
at several values of T is displayed in Figs. || and |(| Fig. || 
shows the thermal conductivity of iron plasma in the den- 
sity range appropriate to the outer envelope of a neutron 
star at lgT[K]= 6.5, 7.0, 7.5 and 8.0. Fig. | shows the 
thermal conductivity of matter composed of H, He, C, Fe 
at T — 10 7 K in about the same density range (related 
to the outer envelope of a neutron star or to the core of a 
white dwarf). We assume that no impurities are present, 
but we include the contribution of ee scattering in addi- 
tion to the ei scattering. The ee scattering contributes 
to v K , i.e, lowers the thermal conductivity. According to 
Paper I, 

_ ee _ 3a 2 (fc B T) 2 /2fcp\ 3 



2tt 3 hm*c 2 V/ctf J 



J(xr,y), 



(32) 



where y = V3T pc /T, T pc = (h/k B ) 
electron plasma temperature, and 



Y / 47re 2 n e /m* is the 



6 

5x 2 
x In ( 1 + 



+ 5x4 ) 
2.81 



y 



3(1 + 0.07414 y) 3 
0.81 v\ 



7T 

IT 



(i3.9i + y y 



(33) 



For a plasma of light elements, Z < 6, the strongest effect 
of the ee collisions takes place at temperatures compara- 



Fig. 5. Thermal conductivity vs density in iron plasma 
at several T. Various symbols show numerical results, and 
curves are the fits. Conductivity due to ee collisions is 
included but is entirely unimportant. Impurities are ne- 
glected. 



ble to degeneracy temperature Tp (Lampe 1968; Urpin & 
Yakovlev 198C), as confirmed by Fig. [| For Fe matter, ee 
collisions are unimportant. For H and He plasmas, their 
effect is more pronounced at lower p, where the chosen 
temperature, T = 10 7 K, is closer to Tp. We conclude that 
ee collisions do not play a leading role in thermal transport 
by degenerate electrons but should be taken into account 
in a plasma composed of light elements. 

The strong effect of chemical composition (Fig. ^ is 
related to the Z-dependence of the collision frequencies 
' K . The higher is Z, the larger is v^ K , and the lower is 
the conductivity. Comparing Figs. || and || we see that a 
temperature variation by a factor of 30 can change the 
thermal conductivity of iron plasma much less than alter- 
ing the chemical composition from H to Fe at fixed T . This 
effect has important consequences for the relationship be- 
tween surface and internal temperatures of neutron stars 
(Paper I). 

Finally let us consider partially ionized matter. 
We expect that this case can be considered in the 
mean-ion approximation (Sect. |3.3| ). Fig. [7] shows the 
dependence of the Coulomb logarithms k a and A K 
on the effective Coulomb plasma parameter T = 
Z 2 ff e 2 (47rn i /3) 1 / 3 (fcB7 1 )~ 1 for a partially ionized Fe mat- 
ter at p = 10 3 g cm~ 3 and Z c g = 4. The assumption that 
Z e ff is independent of temperature is unrealistic, and we 
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Fig. 6. Thermal conductivity vs density in matter of dif- 
ferent chemical composition (H, He, C, Fe) at T = 10 7 
K. Filled circles are numerical values and curves are fits. 
Lower curves and circles are calculated including ei and 
ee collisions, while upper curves and circles take account 
of ei collisions alone. 



adopt it for illustration only. Open and filled symbols show 
numerical values of and A K , respectively, calculated 
in the mean-ion approximation, while dotted and dashed 
lines are the fit curves. Although we did not include the 
data with Z -C A/2 in the fitting, our fit formulae appear 
to be robust and reproduce them quite well. 



5. Conclusions 

We have reconsidered the electrical and thermal conduc- 
tivities of degenerate electrons produced by ei scattering, 
using the modified structure factors of ions as suggested 
in Paper II. We have analysed the electron transport in a 
wide range of densities, from about 10 2 — 10 4 gcm~ 3 to 
- 10 7 - 10 10 gem" 3 , and temperatures T ~ 10 4 — 10 9 K, 
for chemical compositions of astrophysical importance. 
The obtained conductivities differ significantly from those 
calculated previously in a wide range of temperatures, 
Tm/5 Si T ^ 5T m , near the melting temperature T m of 
Coulomb crystals. Our new approach has reduced unre- 
alistically large jumps of the transport coefficients at the 
melting point obtained in the earlier works. This, in turn, 
allowed us to develop a unified description of electron con- 
duction in liquid and crystal matter and obtain an effec- 
tive potential for the ei interaction. 



Fig. 7. Coulomb logarithms versus effective Coulomb 
coupling parameter for a partially ionized iron plasma 
with effective charge Z e g = 4 at p 
mean-ion approximation. 



10 3 g cm 3 in the 



The ei scattering, which we studied in this article, is 
known to be the most important mechanism of electron 
relaxation under prevailing physical conditions in the en- 
velopes of neutron stars and in the cores of white dwarfs. 
We expect that the improved transport coefficients can 
be used to solve various problems of the physics of neu- 
tron stars and white dwarfs (cooling, evolution of accreting 
stars, nuclear burning of matter, pulsation modes, evolu- 
tion of magnetic fields, etc.). 
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Appendix: analytical fits to Coulomb logarithms 

Carrying out the integration in Eq. (ph with the effective 
potential given by Eq. (|25| ) yields the following expressions 
for the Coulomb logarithms, which enter Eq. ((7) for the 
effective collision frequencies: 



A fit = 



Ai(s, to) j A 2 (s,tu) 



(A.l) 
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Here s = q*/(2k F ) 2 , w = w(2k F ) [see Eqs. @-(|l|)], 
s + 1 , 



2Ai («,«;) 
2A 2 (s,w) 



In- 



(l-e"») 



s s + 1 

- (1 + sw) e sw [Ei (sw) - Ei (sw + w)] , (A.2) 



1 + w 



-(1 



2s In 



s + 1 



W 8 + 1' S 

+ s(2 + sw) e sw [Ei (sw) - Ei [sw + w)} , (A.3) 

and Ei(x) — y~ 1 e~ v dy is the exponential integral 
given, for example, by the rational-polynomial approxi- 
mations in Abramowitz & Stegun (1972). 

In the particular case of s <C 1 and s <C w~ x , numerical 
cancellation of large numbers can be avoided with aid of 
the asymptotic expressions 



Ai(fl->0,u;) = -[E 1 (w) + lnw + i\, 



A 2 (s -> Q,u;) 



1 + w 



2w 



where 7 = 0.5772 ... is the Euler's constant. 
In the limiting case of w <C 1, we obtain 



Ai(s,w < 1) 
A 2 (s,w < 1) fa w 



2s + 1 , s 
w I s In — 



1 



2s + 2 
1 - 3s - 6s 2 
4s + 4 



3 . s + 1 
- In 

2 s 



(A.4) 
(A.5) 

(A.6) 
(A.7) 



In the opposite case of w > 1, the familiar expressions 
for the Debye-like screening are recovered: 



1 A s + 1 1 

Ai s,w» 1 = - In — 

2 V s s + 1 



A 2 (s,w > 1) 



2s + 1 s- 

7: ^ — s in — 

2s + 2 



1 



(A.8) 
(A.9) 
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